:arXiv:hep-ph/9901286v3 18 Jul 2000; 



UM-P-98/55, RCHEP98/14 

The Color Singlet Relativistic Correction in J/ip Photoproduction 



C.B. Paranavitane*, B.H.J. McKellar 

Research Center for High Energy Physics, School of Physics, 
University of Melbourne, Parkville, Victoria 3010, Australia 



J.P Mat 

Institute of Theoretical Physics, Academia Sinica, 
P.O. Box 2735, Beijing 100080, China 
(February 1, 2008) 



Abstract 



The 0{y 2 ) relativistic correction for inelastic J/tp photoproduction, in which 
heavy quark pairs are in the dominant Fock state of the quarkonium, is studied 
in the framework of NRQCD factorization. An assessment of its significance, 
particularly in comparison to the color octet contributions, is made. It is 
found that the impact on the energy distribution is negative in certain regions 
of phase space. The predictions are compared with photoproduction data from 
DESY-HERA. 
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I. INTRODUCTION 



Recent theoretical developments have been responsible for a resurgent interest in quarko- 
nia physics: a factorized form for quarkonia production rates , based on the effective field 
theory, Non Relativistic QCD (NRQCD) has been proposed. In this formalism, the 
inclusive production rates can be expressed as a function of the rates for creating pairs of 
free heavy quarks, QQ, in various configurations, and by a series of non-perturbative matrix 
elements. The former may be computed using perturbative QCD, while the matrix elements, 
which describe the transition of the QQ into a hadron, are defined in NRQCD. Contrary to 
previous approaches, NRQCD factorization allows a systematic determination of relativistic 
corrections, which describe the effects of a heavy quark moving at a small velocity v in the 
rest frame of the quarkonium. The velocity scaling rules of NRQCD j0|, which reveal the 
power law behavior with respect to v of the operators of the effective field theory, may be 
used to organize the matrix elements into a hierarchy. Although each quarkonium resonance 
is identified by the QQ having a particular set of quantum numbers (such a QQ state is 
the dominant Fock state of the quarkonium), the relativistic corrections to the production 
rates for that quarkonium, which accommodates the relative motion of the Q and Q, may 
introduce a QQ state with other quantum numbers. For example, the pair may be in a color 
octet configuration. The corresponding channels in the production rate, despite being at 
higher orders in v, can be phenomenologically significant in certain kinematic domains, due 
to the particulars of the production process. 

One of the experiments used to test NRQCD factorization is photoproduction of quarko- 
nia, in which a real photon collides with a nucleon. Charmonium photoproduction has 
been investigated in the framework of the color singlet model ||, and more recently us- 
ing NRQCD factorization ||, whose predictions have been subject to a comparison with 
preliminary data from HERA In addition to the channel in which the QQ is in the 

dominant Fock state of the J/ip, corrections, corresponding to the QQ in spin singlet and 
spin triplet, S and P wave, color octet states, have been considered in the studies. The 
numerical values of the color octet matrix elements, were estimated in an earlier work J?|, 
by fitting theoretical predictions for prompt charmonium production in pp processes, to data 
from the Fermilab-Tevatron. Using these values, however, leads to a significant discrepancy 
with the photoproduction rates at HERA: the contributions from octet channels overwhelm 
the data in certain regions of phase space. 

A number of uncertainties which enter into calculations may be responsible for the incon- 
sistency between the HERA and Tevatron data for charmonium production, when viewed 
from the perspective of NRQCD. Among these are uncertainties in parton distribution func- 
tions, the value of the charm quark mass, higher twist effects [II], and higher order pertur- 



bative QCD corrections. In the work by Kniehl ||20|| , important higher order perturbative 
QCD effects, due to initial state multiple gluon radiation, were included when analyzing 
hadroproduction at the Tevatron. The re-estimated matrix elements have values signifi- 
cantly smaller than the originally determined ones. Correspondingly, the HERA data is in 
much better agreement with the NRQCD predictions. This is true throughout most of phase 
space, except near certain boundaries where the octet contributions diverge. The authors 
in fT2fl, have suggested this is due to the non- relativistic expansion breaking down in such 



regions. They have proposed a resumation of the non-relativistic expansion, which requires 
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the introduction of universal structure functions. 

The contributions to J /if) photoproduction from the color octet channels are at order 
v 4 in the small velocity expansion, relative to the leading order, color singlet term. There 
is a relativistic correction at order v 2 , incorporating the effects of the relative motion of 
the heavy quark pair in the dominant Fock state of J /if). The purpose of this paper is 
to examine this correction in the framework of NRQCD factorization, which requires the 
introduction of an additional matrix element. Previously, models have been constructed 
to describe such an order v 2 effect. The one parameter model presented in |L8| has been 



used to examine quarkonia photoproduction by Jung et al. [O; a two parameter model has 



also been proposed |14]]. The 0(v 2 ) contribution has not been computed within NRQCD. It 
has been pointed out by Bodwin et al. ||, that corrections of 0{y 2 ) may be as important 
as radiative corrections of 0(a s ), which have already been calculated for photoproduction 
T0fl , an d found to be significant. Without a detailed understanding of the 0(v 2 ) term, it is 



not possible to make an assessment of the significance of the color octet terms: the 0(v 2 ) 
contribution is clearly necessary for consistency, according to the power counting rules of 
NRQCD. 

The remainder of the paper is organized as follows. In the next section, some details of 
the calculation of the relativistic correction is described. In section [TTT| , numerical results, 
indicating the impact of the v 2 correction, are presented, together with a comparison with 
data from HERA. Section fV] contains some concluding remarks. 



II. PHOTOPRODUCTION 

NRQCD factorization allows the inclusive cross section for the photoproduction of J /if) 
to be written in the following form: 

da_ 1+i ^j/,p +x \ - F n /snJ/iP\ m 

di iW 1 - 4 

The coefficients, F n depend on the kinematic invariants of the process in which a photon 
7 and a parton i from the initial hadron, react to produce a QQ. Here m is the mass of 
heavy quark. The production is a short distance process, and is therefore calculable using 
QCD perturbation theory. The subscript n = 25+1 L^j' 8 ^ specifies the spin S, orbital angular 
momentum L, and the total angular momentum J, of the QQ. The 1 and 8 means that 
the pair may be in a color singlet or color octet state. The matrix elements (O/J^), of 
mass dimension d n describes the transition of the QQ in the state n, into a J /if). Although 
J /if) is identified as a QQ mainly in a 3 Si configuration, the short distance process may 
produce a heavy quark pair which is not in a 3 S[^ state. The formation of the bound 
state can involve the emission of soft partons, which return the pair into the appropriate 
spectroscopic configuration. The relative importance of various terms in the sum in Eq. (|IJ), 
may be determined by how each term is suppressed by powers of the strong coupling constant, 
a s , whose dependence lies in F n , and the velocity of the heavy quark in the rest frame of 
the quarkonium v, appearing in (O/J^). For J /if), v 2 ~ 0.3. 

The dominant short distance contribution to inelastic photoproduction of J /if) is the 
photon-gluon fusion process 
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1 (k)+G(g 1 )^QQ n (P) + G(g 2 



(2) 



The four momenta are specified in parenthesis. The process (§) at leading order in QCD 
perturbation theory is represented in Figs. 1(a) and 1(b). Both diagrams allow the heavy 
quark pair to be produced in a color octet state, while the QQ can only to be in a color singlet 
state in Fig. 1(a). Subprocesses involving light quarks in the initial state are considerably 
suppressed at the energies available at HERA 0. Resolved processes, in which the photon 
breaks up into partons, which subsequently interacts with the nucleon, have been analyzed 
within NRQCD factorization [HJ. These are important in the kinematic regime z < 0.3, 



where z = Ej/^/E 1 , in which the energies are defined in the nucleons rest frame. 

The formation of J /if) from a QQ at leading order in the small velocity expansion is 
described by the matrix element 



J/*(3 S (l> 











(3) 



The field if) annihilates a heavy quark, while x creates a heavy anti-quark. The matrix 
element is proportional to the rate of production of a J/if> from the NRQCD vacuum, via 
an intermediate state in which the heavy quark pair is in a 3 S[^ configuration. S represents 
light hadrons in the final state, and J z is the helicity of the J /if), both of which are summed 
over. 

The relativistic correction, appearing at 0(v 2 ) relative to the leading contribution, in 
which the heavy quark pair is also in the dominant Fock state of the J /if), involves the 
matrix element 

2 



(1)- 



pJ/V>(3£( 



where x 1 D ip 
the form, 







2 D 



X + h.c 







(4) 



X^(D)-?/> — (Dxyif)- Therefore, the differential cross section to 0(v 2 ) takes 



da 
di 



f( 3 ^ 1} ; 



G( 3 s[ iy ) 



rrr 



(5) 



Of the octet contributions which appear at 0(v A ), the 1 si S) and 3 Pj 8) channels were found 
to be the most important. The structure of the matrix elements which describe these 
transitions, are similar to Eqs.flU) and (|), except for the operators between the quark 
and anti-quark fields; they may be found in the literature 0. The aim in this section is 
to determine the coefficient G( 3 5'j 1 ' ) ), containing the short distance information associated 
with the v 2 correction. 

The low energy equivalence of full QCD and NRQCD allow the short distance coefficients 
F n to be determined by the following matching condition: 



da 
di 



1 +G^QQ n +X 



E 



pQCD 



dn-4 



m 



(0\O^\0) 



pNRQCD 



(6) 



On the left hand side of Eq. |6|), the cross section for the production of a free QQ is computed 
using perturbative QCD. The matrix elements on the right hand side are evaluated using 
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perturbative NRQCD. Both sides are expanded in powers of q, the momentum of the heavy 
quark in the rest frame of the quarkonium. The F n are then found by comparing powers 
of q. The QCD side of the matching condition involves expanding the amplitude for the 
process, Eq. (H), in powers of q. The leading terms contain QQ's in the dominant Fock 
state of the J/ip, 3 5i , as well as some octet configurations: 3 S[ 8 ^ and 1 S^\ Terms at 
order q in the amplitude project out heavy quark pairs in 3 Pj 8 ' ) states, but do not contain 
a heavy quark pair in 3 S[ 1 ' > states. In order to determine the 0(v 2 ) contribution in (j^), it is 
necessary to expand the QCD amplitude to 0(q 2 ): at this order, heavy quark pairs in 
configurations are present. The reason that the color octet contributions appear at a higher 
order in v, in (H), despite being at a lower order in the perturbative amplitude, is because 
the power counting is different for hadronic states: additional v suppression appears in the 
hadronic projection operators of the color octet matrix elements. 



Some conventions used in []T6| will be adopted in this work. The four momenta of the 
heavy quark and anti-quark is 

P" = ^~ A 1V, (7) 
where A is the Lorenz boost matrix from the rest frame of the J/i/j, to the frame in which 



it is moving with four momentum P = (y4(m 2 + |q| 2 ) + P 2 , P). The partonic Mandelstam 

invariants are defined in the standard manner: s = (k + g{) 2 = (P + g 2 ) 2 , t — (k — P) 2 = 
(gi — g 2 ) 2 and u = (g\ — P) 2 = (k — g 2 ) 2 '. Only two of these invariants are independent, their 
relationship is given by 

s + i + u = P 2 = A{m 2 + |q| 2 ). (8) 

We choose u to be expressed in terms of the other two invariants, and |q| 2 , which is associated 
with (V J ^( 3 S[ 1 " > )). Two contraction identities found in reference [16] will be frequently used 
in the calculation: 

pfj, pu 

Af Ar = -g, v + -j^. (9) 

First we will consider the QCD side of Eq. The transition matrix for the process 
may be written, 

/ ly+G-*QQ n +G = e,(k)sl( gi )sf(g 2 )u{p)T^v{p). 
The color indices of the gluons are specified by a and b. Following the techniques used 



m 



HI) i s decomposed in terms of the sixteen basis matrices 1,75,7^,7^75, and cx 



When these are inserted between the quark and anti-quark spinors u(p) and v(p), T may 
be expressed in terms of the color triplet, Pauli spinors, £ and rj. Projecting out the color 
singlet channel and expanding to 0(|q| 3 ), the transition matrix takes the form 



6 



r = e»(k)el( gi )ef(g 2 )^(ar P + C^Y)^'^ + ... (10) 

The terms af Q/3 and bf^ are functions of the kinematic variables and the Lorenz boost 
matrices. The ellipsis represents contributions in which the QQ pairs are in color singlet P 
wave states, spin singlet states, or terms which contribute to processes occurring at 0(v 3 ) or 
higher. To obtain the rate, T is multiplied by it's complex conjugate, and the appropriate 
sums and averages over spin and color degrees of freedom are made. Before the momentum 
space integration is carried out, a change of variables is performed: (p,p) — > (P,q), to 
facilitate the matching as defined in Eq. (|6]). The perturbative cross section becomes 



dd_ y+G -, QQn+x L_ ( _L ^ r rf3 q 



dt 16tts 2 V2 2 .4iV c 2 J J (2?r) 



A--+ fh {2) - | n | 2 4- r (2) n m n n 



a 1 ®a 3 + ... 



(11) 



For convenience, the notation suppressing the Pauli spinors and the sums over the spins has 
been employed 0. The coefficients appearing in the transition matrix are related to the 
ones appearing in (|TT|) . Using the identities in the contraction of the Lorenz indices 
relates the A^, bff and c\f mn , to af a/3 and b^ n . 

or /, a^jj = Ay + 6g ) |q| 2 + 0(|q| 8 ) 

fj,a/3,imn ' u imn u ' fia/3,j ^ijmn ' ^VIVll/ 

It must be kept in mind that the factor A{j implicitly retains some q dependence which 
will become fully exposed once contractions over the remaining indices are carried out. The 
term Aij/2m 2 is a consequence of the change of variables. Using rotational symmetry, Eq. 
( p|) takes the form 



da 



- < :>■>„ v ... >_(_!. \ t JLJL i, M, 7 r < .r n J . i ( 3c (2) + 6 (2) _ £*L_ I ,12; 



dt 16tts 2 \2 2 .4iV 2 / 7 (2tt) 3 3 m 3 I " 2m 2 




Only the leading order terms in |q| 2 have been kept in the contractions of c^ mn and bf^ 1 . 
q° and cj 2 ^ are defined to be 

A.. _ _(0) , r ( 2 )| n |2 

Equation ([TJ) displays all the q 2 dependence necessary to determine G( 3 Si). The coeffi- 
cients have been expressed in terms of the invariants s and t using the contraction identities. 
They are somewhat lengthy and will not be presented here. 

Next the NRQCD side of the matching condition is considered. The QQ differential cross 
section, only including the contributions we are interested in, may be written 

E^l {0 ^( 3 si 1] )) + ^p-iv^es^)} + 

m 2 m 4 
ci/3 c( 1 ) 3 r,( 1 )\ _ 

m 4 h v^fsP,' D?)) + ■■■ (13) 
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It is necessary to take into account another color singlet matrix element with d n = 8, 



3 s?\ 3 d?] 





1 


<• 






2 



. (14) 



The symmetric traceless tensor, formed from the covariant derivatives is specified by D D = 

(dT) + S 5 )/2- D • D /3. Here (V QQ ( 3 Si, 3 D x )) describes the overlap of the transition 
amplitudes in which the QQ is in spin triplet S and D wave states. Therefore, F( 3 Si 3 Di) 
mixes with F( 3 S\) in perturbation theory. (V®®( 3 Si 3 -Di)) appears at 0{v 2 ) in perturbation 
theory, where quarks and gluons are free. However, because the velocity counting is different 
for hadronic states, the corresponding hadronic matrix element (V J ^( 3 Si 3 -Di)) is at order 
v 4 . This is because it requires two chromo-electric dipole transitions to take a QQ in a 
3 S[^ configuration, to a 3 D^ state, costing a factor of v 2 in the projection operator of 
(V J ^( 3 S[ 3 )}. The covariant derivatives in the operator contributes another v 2 . 

Using a non-relativistic normalization convention, it is possible to realize the right hand 
side of (El) at leading order in perturbative NRQCD: 



d 3 q 



3q(1)> 



F( 3 S\ 



1 



(2tt) 3 
+ 3 \G( 3 S{ 1] 



-a 



m z 



3m 4 



3F{ 3 S{ 1)3 DP 



Iq.cr ® q.cr 



F( 3 ^ (1) 3 



'1 5 



0-h 



qV 



(15) 



One may worry about the use of different normalization of states on either side of (0). 
However this is not a problem as only physical quantities are being considered, and the 
Pauli spinors, which are not shown explicitly, have the same normalizations: £t£ = 1 and 
ryr] = 1. Braaten et al. [[UJ defined the matrix elements with relativistic normalizations, 
which were subsequently related to to the matrix elements defined in NRQCD. 

By comparing Eq. ( |I~5"D with Eq. flT2|), we finally obtain the short distance coefficients. 
The leading order color singlet contribution is well known: 



F( 3 S[ 



1 32m 
16ns 2 27 



(47ra s ) 2 (47rct)e 2 



s 2 (s - Am 2 ) 2 + i 2 (i - Am 2 ) 2 + (4m 2 - s - i) 2 (s + if 
(s-4m 2 ) 2 (t-4m 2 ) 2 (s + t) 2 



(16) 



The strong coupling constant and the fine structure constant are given by a s and a. The 
electromagnetic charge of the charm quark in units of e, is e c . As a consequence of Eq. (Bj), 



Eq. 



is expressed in terms of the combination 4m 2 — s — i, instead of u. The perturbative 



coefficient of the relativistic correction for the color singlet channel is found to be 

G( 3 S{ 1] ) 



32(4vra) 2 (47ra)e 2 (-6144m lu s 2 + 1280m 8 s 3 + 960m 6 s 4 



rrr 



167TS 2 

- 336m 4 s 5 + 28m 2 s 6 



4096m 10 st - 512m 8 s z t + 3200m b s 3 t 



8s2+ 



,6 -3i 



- 1008m 4 sH + 72m z sH + s n t - QU4m w t z - 512m 8 st z + 3200m°sH 



,2-5 



10+2 



,6-2+2 



- 1408m 4 s 3 t 2 + 92m z sH z + 3s°t z + 1280m 8 t d + 3200m°st 



,2-4+2 



5+2 



,8+3 



,6;+3 
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- 1408m 4 s 2 F + 112m 2 s 3 i 3 + 5s 4 P + 960m 6 t 4 - 1008m 4 st 4 

+ 92m 2 s 2 t 4 + 5s 3 t 4 - 336m 4 t 5 + 72m 2 si 5 + 3sH 5 + 28m 2 t 6 + si 6 )/ 



(81m(s - 4m 2 ) 3 (t - 4m 2 ) 3 (s + *) 3 ) 



2\3/ 



(17) 



This result is different to the one found in [13|]. First, the relativistic correction in this 



work is expressed in terms of a matrix element defined in NRQCD; in jl3|, it is correspond- 
ingly a parameter e, defined via 



M 



J/tP 



2m + e, 



(18) 



where e was taken to be positive. As pointed out in [|TTJ, this parameter can be negative 
and can be defined as a matrix element in NRQCD which is different to the one appearing 
in Eq. (5). Secondly, the correction calculated in Hl3| , is based on the model proposed in 
|T8|| , in which the contribution of the momentum q to the energy of the heavy quark is 
neglected-that is the energy is taken to be Mj/^/2, in the quarkonium rest frame; it really 
should be y/ m 2 + q 2 . This implies that our result cannot be obtained from that presented 
in [13| by setting Mj/^ = 2m, and the parameter e proportional to (V H ^'( 3 S[ lS) )}. 



III. NUMERICAL RESULTS 

In this section an assessment of the significance of the order v 2 correction shall be made, 
using the analytic result (|17|). The predictions will be compared to photoproduction data 
from HERA. In ||, the color octet channels 1 Sq 8 " > and 3 -Pj 8 ' ) were found to be important in 
photoproduction of J/ip. The corresponding matrix elements were taken to have numerical 
values consistent with the linear combination [0] 

/QWPsP)) + ^(£> J/ ^( 3 P (8) )> = 6.6 x 1(T 2 GeV 3 , (19) 
m 2 

determined in the original analysis of hadroproduction rates at the Tevatron. However, these 
values lead to predictions which are significantly greater than the experimental results. The 
inclusion of effects due to initial state gluon radiation modifies the linear combination of the 



matrix elements thus [|20 



(O j ^( x 4 8) )> + ^(O J/V, ( 3 P (8) )> = 5.27 x 1CT 3 GeV 3 . (20) 



In this study, we shall assume the latter, and take (O^^S^)) = (C J /^( 3 P (8) ))/m 2 = 
1.2 x 10~ 3 GeV 3 . The leading order color singlet matrix element is taken to have the value: 

(O j/ ^( 3 S{ 1] )} = 1.3 GeV 3 . (21) 

A naive estimate for (V J ^( 3 S[^)} shall be used, based on the power law suppression of 
(pJ/tf(3 5 (i))) re i at i ve to (O j ^( 3 S{ 1] )}: 

(V^^CS^)) = 0.85 GeV 5 . (22) 
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The other inputs for the computations are as follows: the strong coupling constant at the 
charm mass scale is taken to be a s (m) = 0.3, with m = 1.48 GeV. The MRSA determination 

of the gluon structure function is used JT7], with Aq J ^ /_4 = 300 MeV, and Q 2 = (2m) 2 . 

The total cross section as a function of the center of mass energy is considered first. The 
constraints: z < 0.8 and pt > 1 GeV, valid at HERA, are adopted. The K factor, consistent 



with these cuts, from the NLO corrections to the color singlet channel |10] is K ~ 1.8. 
Figure 2 shows the total cross section as a function of the photon-nucleon center of mass 
energy, yfs. The leading order color singlet contribution (long dashed line), accounts quite 
well for most of the data from HI || (circles) and ZEUS |4| (squares). The v 2 correction 
(short dashed line), together with the color octet terms (dot-dash line), contributes only 
marginally. The octet contributions are greater than the 0(v 2 ) correction by approximately 
a factor of four, for large center of mass energies, with the values of the matrix elements 
that have been adopted. If the octet matrix elements were chosen to be consistent with 
(P3|), then the v 2 correction would be smaller by approximately a factor of ten, relative to 
the octet contributions. However, (O j ^( 3 Si)) may be somewhat larger than 0.85 GeV 5 , 
therefore both corrections may have similar magnitudes. 

Next, the da/dz distribution is considered, as a function of z. Figure 3(a) shows the 
ratio of the v 2 to the octet contributions, at \fs = 100 GeV (solid line) and \/s = 14.7 GeV 
(dashed line). At y/s = 100 GeV the cut, p t > 1 GeV, consistent with results from HERA 
is adopted, while at y/s = 14.7, p t > 0.1 GeV, valid for the European Muon Collaboration 
(EMC). For \/s = 100 GeV, it can be seen that the the v 2 channel dominates at low z. 
For z > 0.63 the ratio is negative, indicating that the v 2 distribution is negative in that 
region. At z ~ 0.7 the ratio has a minimum value of about 0.45, meaning that in this 
region, the v 2 contribution is significant with respect to the color octet channels. Again, if 
the original determinations of the octet matrix elements were adopted, the v 2 contributions 
will give a negligible result at intermediate and large values of z, while at low z, they would 
be comparable to the color octet contributions. At the EMC energy regime, \fs = 14.7 
GeV, v 2 corrections are relatively more substantial at low z, while at larger values, not as 
important. Figure 3(b) is the z distribution at the center of mass energy y/s = 100 GeV. 
The long dashed line represents the leading order contribution. The sum of the leading 
order and v 2 correction are shown with a dot-dash; the sum of the leading order and octet 
contributions is shown with the dashed line. The net result is shown with a solid line. The 
circles represent data from the HI Collaboration ||. The v 2 correction tends to decrease the 
distribution for larger values of z, while at smaller values, it is enhanced. Figure 3(c) shows 
the distribution at y/s = 14.7 GeV, displaying a similar behavior, although the v 2 effect is 
not as pronounced. The data is from the EMC ||. It should be noted that the large z region 
where the v 2 contribution is most substantial, is also where the non relativistic expansion 



breaks down |12]. Our analysis, leads to a qualitatively different behavior for the energy 
distribution, than in the work by Jung et al. ||13|| . In their case there is an enhancement of 
da/dz for large values of z. 

The transverse momentum distribution for photoproduction was also examined. How- 
ever, it was found that the v 2 correction leads to a negligible contribution to the overall 
differential rate. It is expected that a similar v 2 correction appears in p + p —>■ J/ip + X 
processes at the Tevatron, where color octet fragmentation channels dominate at large trans- 
verse momentum. The v 2 correction can be obtained by realizing that it originates via a 
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process similar to that depicted in figure 1(a), except with the initial photon replaced by 
a gluon. The color singlet rates, including the v 2 correction is obtained by multiplying 
F( 3 S^) and G( 3 S[ ) by the factor 5a s /(96ae^), which takes into account the differences in 
the initial state. However it was found that the v 2 correction contributes insignificantly to 
the transverse momentum distribution for moderate transverse momenta, at the Tevatron 
center of mass energy y/s = 1.8 TeV. 



IV. CONCLUSION 

In this work, the relativistic correction to the color singlet channel in J / ip photoproduc- 
tion was calculated using the NRQCD factorization formalism. The correction is described 
by an additional matrix element which appears at 0(v 2 ) relative to the leading order color 
singlet term, according to the velocity scaling rules of NRQCD. This work contrasts with 
previous studies in that it draws qualitatively different conclusions about the behavior of 
observable. Specifically, the contribution leads to a decrease in the energy distribution for 
large values of z, while for smaller values there is an enhancement of the differential cross 
section. The effect of the correction in photoproduction seems to be more significant for 
larger center of mass energies. In some regions of phase space, the v 2 correction is com- 
parable to the color octet contributions, with the values of the parameters that were used. 
However, it has been noted earlier that the rates are sensitive to a s (m), m, and the choice 
of parton distribution functions, which are not precisely known. Certainly it seems that the 
uncertainty in these parameters, correspond to variations in the rates which are of similar 
magnitude, or greater than the effect of the 0(v 2 ) correction [10|. Therefore, at the present 
level of precision, it is not possible to observe the effect of the correction in experiments. 
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Figure Captions 

Figure 1: (a) One of six Feynman diagrams for the production of a QQ pair. The others 
are obtained through permutations of photon and gluon lines. The QQ may be in a color 
singlet or color octet configuration, (b) A typical diagram in which the QQ's are only in 
color octet states. 

Figure 2: The total cross section a{\fs) (nb) as a function of the center of mass energy 
y/s (GeV). The leading order color singlet contribution is shown with a long dashed line. 
The 0{y 2 ) correction is the dot-dash line; the octet contribution is the short dash line. The 
total rate is the solid line. 

Figure 3: (a) da/dz: the ratio of the v 2 contribution to the octet contributions as a 
function of z. The solid line is at ^/s = 100 GeV, and the dashed line is at ^/s = 14.7 GeV. 
(b) da/dz (nb) as a function of z, at the center of mass energy y/s = 100 GeV (HERA). The 
long dashed line is the leading order color singlet contribution. The sum of the leading order 
and 0{y 2 ) channels is the dot-dash line. The sum of the leading order and octet contribution 
is the dashed line. The solid line represents the total rate. The circles are the data from the 
HI collaboration, (c) da/dz (nb) at ^/s = 14.7 GeV. The lines represents contributions as 
described in (b). Data from EMC are shown with solid circles. 
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